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Abstract 

Studies on the genetic diversity and relatedness of zoo populations are crucial for Implementing successful breeding 
programmes. The European wildcat, Felis s. silvestris, is subject to intensive conservation measures, including captive 
breeding and reintroduction. We here present the first systematic genetic analysis of the captive population of Felis s. 
silvestris in comparison with a natural wild population. We used microsatellites and mtDNA sequencing to assess genetic 
diversity, structure and integrity of the ex situ population. Our results show that the ex situ population of the European 
wildcat is highly structured and that it has a higher genetic diversity than the studied wild population. Some genetic clusters 
matched the breeding lines of certain zoos or groups of zoos that often exchanged individuals. Two mitochondrial 
haplotype groups were detected in the in situ populations, one of which was closely related to the most common 
haplotype found in domestic cats, suggesting past introgression in the wild. Although native haplotypes were also found in 
the captive population, the majority (68%) of captive individuals shared a common mtDNA haplotype with the domestic cat 
{Fells s. catus). Only six captive individuals (7.7%) were assigned as wildcats in the STRUCTURE analysis (at K = 2), two of 
which had domestic cat mtDNA haplotypes and only two captive individuals were assigned as purebred wildcats by 
NewHybrids. These results suggest that the high genetic diversity of the captive population has been caused by admixture 
with domestic cats. Therefore, the captive population cannot be recommended for further breeding and reintroduction. 
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Introduction 

Due to the ongoing threats to biodiversity, species conservation 
is still a challenging task [1]. Although in situ conservation 
represents by far the most efiective way to protect endangered 
species, it is evident that not all species can be preserved in their 
natural habitats. Therefore, ex situ conservation and reintroduc- 
tions have become common measures of species conservation 
[2,3]. Despite the great advances in organization and logistics 
(studbooks, online databases) and increasing expertise in veterinary 
medicine, inbreeding and outbreeding stUl pose serious threats for 
ex situ conservation [4-7] . Captive breeding may be compromised 
by unknown founder origin and relationships [8], undetected 
hybridisation [9], limited representation of natural genetic 
diversity [10] or studbook errors [11,12]. These problems can 
affect both the genetic integrity of the captive population as well as 
the success of reintroduction projects. Therefore, data on the 
genetic diversity and relatedness in zoo populations can provide 
valuable information for improving conservation programmes 
[13]. 

The influence of captive breeding on the genetic diversity of 
endangered species has already been subject to a number of 
theoretic and genetic studies [7,14-16,17]. However, the majority 



of the genetic studies focus on highly endangered species, for 
which studbooks or international breeding programmes already 
exist [18-20]. This circumstance has two drawbacks: (1) Species 
with small population size might have a reduced genetic diversity 
even in their wild populations. Therefore, the genetic diversity of 
captive populations of such highly endangered species might not 
reflect the effects of captive breeding. (2) Studbooks are usually 
kept in a way to avoid the loss of genetic diversity or inbreeding 
depression. Hence, it is likely that genetic analyses of ex situ 
populations will confirm the success of such coordinated breeding 
schemes [16,21]. However, studbooks exist only for ca. 850 species 
(WAZA 2010), while ca. 10.000 species are currently kept in zoos 
and bred without any coordination (ISIS 2010). Many of these 
species are still rather common in the wild, but nevertheless some 
of them are regionally endangered and sometimes even subject to 
reintroductions. In order to assess the influence of uncoordinated 
breeding on the genetic diversity of a captive population, it is 
useful to study one of these species rather than well managed 
captive populations. Comparing the genetic variability of ex situ 
populations and wild populations also requires that a target species 
has viable wild populations and a sufficiendy large captive 
population. Furthermore, it should have been successfully bred 
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for many generations in captivity as genetic erosion may 
accumulate over time [18]. 

The European wildcat {Felis silvestris silvestris) is an ideal study 
object for such an analysis. It is still rather common in Europe and 
classified as Least Concern on the lUCN Red List of Threatened 
Species [22]. However, the populations of this species are 
fragmented and in some countries declining [22]. Wildcats are 
kept and bred in many European zoos and captive individuals 
have already been reintroduced in three regions in Germany [23]. 
However, to date no coordinated breeding programme or 
studbook exists for the European wildcat. Hence, there is no 
information available on the captive stock concerning founder size, 
founder origin, captive population size and relatedness. As it is 
planned to establish a studbook for the European wildcat in the 
near future (A. Sliwa pers. comm.), a genetic analysis of the 
breeding stock can provide important basic information for this 
studbook. Another challenge in the conservation of the European 
wildcat is hybridization with domestic cats, which is a serious 
threat in some wild populations [24-26]. Thus, an additional 
advantage of a genetic analysis prior to coordinated breeding is the 
exclusion of potential hybrids. 

Here, we present the first systematic genetic analysis of the 
European captive population of Felis s. silvestris. We used both 
microsateUites and mitochondrial DNA (mtDNA) sequencing in 
order to assess the genetic variability, structure and integrity of the 
ex silii population. We also analysed samples from a wild 
population to evaluate if the genetic variability in the captive 
population is comparable to the diversity found in the wild. 

Methods 

Sampling and DNA extraction 

We contacted 124 zoos, 64 of which confirmed keeping 
European wildcats. We obtained 80 samples (12 buccal swabs, 8 
tissue samples from deceased individuals, 10 blood samples and 50 
hair samples) from ,30 zoos spanning six European countries (see 
Text S2 in File SI and Tab. S4 in File SI for an overview). All ex 
situ samples were taken during routine veterinary treatments (like 
vaccination, or the placement of transponders for juveniles) or 
whenever a veterinary treatment was necessary. The collection of 
hair samples does not pose a severe stress or hurt to the cats, the 
same is true for buccal swaps from young individuals. Blood 
samples were only takc-n from anaesthetized individuals. 

The in situ sample consisted of 89 individuals (26 tissue samples 
and 63 hair samples) fi-om a natural wild population in the Harz 
mountains (Germany). The hair samples were collected in the 
context of a radio telemetry study (M. Gotz). The tissue samples 
were from a monitoring program for roadside casualties. Both 
studies were carried out with the permit of the Landesamt fiir 
Umweltschutz Sachsen-Anhalt which included the collection of 
samples from European Wildcats (see Text SI in File SI for 
further information). Additionally, we received six samples (3 tissue 
and 3 hair samples) from road kills and live traps collected in the 
wild populations in Rhineland-Palatinate and Saarland. 

Furthermore, we sampled 33 domestic cats with the permission 
of their owners, in order to detect hybrids. 15 of the domestic 
samples came from a private animal shelter which was located 
near the in situ population. We chose this shelter, as it would 
provide a realistic insight into the potential influence local feral 
cats might have on the wild population. For the remaining samples 
we contacted private owners. The private owners tore out a small 
bunch of hair, so that the sample included the follicles. This bunch 
was then sent or given to us in an envelope. 



DNA was extracted from buccal swabs, tissue and blood using 
the DNeasy blood and tissue kit (Qiagen) following the manufac- 
turer's instructions (with special protocols for blood and buccal 
swabs as provided by the manufacturer). For the hair samples, we 
extracted DNA with a modified Chelex 100 protocol, using a IQ% 
Chelex solution with addition of 7 \A Proteinase K (18 mg/ml) per 
sample and overnight lysis [27,28]. No DNA could be extracted 
from seven of the 208 samples (5 from wild populations and 2 from 
the captive population). The respective samples from the wild 
population were hair samples which had been taken several years 
ago and had then been stored at room temperature. The captive 
samples had probably been conserved in denatured ethanol. 

Genotyping 

The six samples from the wild caught individuals in Saarland 
and Rhineland-Palatinate were not genotyped, as the samples size 
from these populations was too small for population genetic 
analysis. They served as reference samples for the mtDNA analysis 
(see below). Therefore only 194 of the 201 samples were genotyped 
at ten di-nucleotide repeat microsatellite loci (FCA08, FCA23; 
FCA43; FCA58, FCA77, FCA78, FCA90, FCA126, FCA132 and 
FCA149) characterized in the domestic cat [29,30]. Amplification 
was performed in a Mastercycler (Eppendorf) using the 
2.5 X5PRIME HotMasterMix (5PRIME). For each PGR we used 
5 |il reaction mix containing 1.2 |xl genomic DNA, 2.2 |J.l 
HotMasterMix, 2.2 |tl water and 0.1 nl forward and reverse 
primers. The PGR conditions were as recommended by the 
manufacturer, with an annealing temperature of 55°G for most 
primers (exceptions: FGA90: 60°C; FCA78: 50°C and addition of 
2.5 mM Mg^""). The 5'-end of each forward primer was labelled 
with a fluorescent dye, either 5-FAM, TAMR.\ or JOE. The 
products were genotyped on a MegaBAGE 1000 automated DNA 
sequencer (GE Healthcare). Fragment lengths were determined 
using Fragment Profiler 1.2 (Amersham Biosciences). To minimize 
genotyping errors due to low DNA concentrations (especially in 
the hair samples), applied a multiple tube approach as 

recommended by Taberlet et al. [31]. Samples were only included 
into further analyses, if they yielded unambiguous results in three 
independent replications (190 samples met this criterion. Tab. S4 
in File SI). Two captive samples had identical genotypes in all 
three replications and it turned out that this individual had been 
transferred to another zoo and sampled twice. Therefore, the 
duplicate sample was excluded from further statistical analysis and 
this left 189 samples for analysis. 

Sequencing 

For 70 wild, 66 captive and 33 domestic individuals we 
sequenced the mitochondrial DNA fragment NADH dehydroge- 
nase subunit 5 (ND5) based on the study by DriscoU et al. [32]. 
Most juveniles of sampled mothers and siblings were excluded and 
later assigned to the same mtDNA haplotype as all sequenced 
juveniles had the same haplotype as their dam. We finally obtained 
88 sequences of wild, 77 of captive and 33 of domestic individuals 
(Tab. S4 in FUe SI). The primers GD-ND56-F1G and GD-ND56- 
R4 [32] were used for amplification in a Multigene Gradient 
Thermal Gycler (Labnet). We used the 5PRIME HotMasterMix 
(5PRIME) for amplification. The PGR product was purified using 
the High pure PGR product purification kit (Roche) according to 
the manufacturer's protocol. Sequencing was performed with the 
DYEnamic ET terminator cycle sequencing kit (GE Healthcare) 
for sequencing reactions run on a MEGABAGE 1000 automated 
sequencer (GE Healthcare). Base-calling was performed in 
Sequence Analyzer 4.0 (Amersham Biosciences). 
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Data analysis 

Genotyping. As the integration of hybrids and genetic sub- 
groupings were unknown, we first searched for genetic structure in 
our samples using Structure 2.3.4 [33]. We assumed admixture 
between groups and used the correlated allele frequency model 
with a burn-in period of 100,000 simulations, followed by one 
million Markov chain Monte Carlo simulations. Tests were run for 
K=l-15 with ten iterations for each K. In order to detect 
potential hybrids between wildcats and domestic cats, we used the 
STRUCTURE results at K = 2 and assigned all individuals either 
to the cluster "wildcats", "domestic cats" or "hybrid" [34]. The 
"hybrid" cluster was chosen, if Q_ values varied between 0.2 and 
0.8 as recommended by Randi and Schulte et al. [35,36]. The 
optimal values for K were assessed using both the method 
described by Pritchard et al. [33], and the method suggested by 
Evanno et al. [37]. The method described by E\'anno et al. [37] 
tends to result in low K values [38,39] and generally works better 
for scenarios with strong genetic differentiation [40] as it detects 
the highest level of differentiation. The method described by 
Pritchard et al. [33] might lead to inflated K-values if closely 
related individuals exist in the sample. Therefore, the clusters 
might more likely represent family groups and lineages [41]. In 
our case the detection of family groups and breeding lines is a 
desired outcome, which provides deeper insight into the genetic 
structure of the captive population. We, therefore, present the 
results obtained by the method described in Pritchard et al. [33]. 
However, we stopped increasing K when Q, values for the next 
cluster dropped below 0.9 in all individuals as proposed by Schulte 
et al. [35]. The individuals were assigned to genetic clusters using 
the highest assignment probability. The dataset was checked for 
nuU alleles using Micro-Checker 2.2.3 [42]. We checked all loci 
for potential linkage using the linkage maps provided by Menotti- 
Raymond et al. [30,43-45] . A test for linkage disequilibrium was 
performed in Fstat 2.9.3.2 using a log-likelihood ratio G-statistic 
with Bonferroni corrections [46]. For the latter test, we removed 
juveniles from the dataset, as these might erroneously suggest the 
presence of a linkage disequilibrium due to their close relationships 
[47]. 

Measures of genetic differentiation have lately been subject to a 
broad discussion [48,49]. Thus, we estimated genetic differentia- 
tion between the genetic clusters with three different estimates, 
Dest, Est and RgT using GenAlEx 6.501 [50,51]. F^t is stiU a 
useful measure if the population split is rather recent (e.g. for 
comparing differentiation among captive populations), while Rgx 
is more useful if the split is deep and if it is likely that mutations 
have already contributed to population differentiation (e.g. for 
comparing wildcats and domestic cats) [52,53]. When we tested 
Dest estimates for correlation with Rgx and Egx, we found a strong 
positive correlation between D^st and Est (R^ = 0.47, P = 0.004) 
but a negative correlation with Rst (R^ = 0.41, P = 0.01). 
Therefore, only the results of Esi and Rs i are provided. Eor 
the analyses of population differentiation, excluded all 

individuals which could not clearly be assigned to a cluster at 
K = 2 based on qi>0.8 and where thus assigned as hybrids (see 
above; all excluded individuals are marked with * in Table S5 in 
File SI). Nested AMOVAs based on Fsx and Rsx were performed 
in GenAlEx with the sample categories "captive", "domestic" and 
"wild" as regions and the six genetic clusters as populations. 

To analyse the degree of hybridization between domestic cats 
and wildcats, we performed an analysis in NewHybrids 1 . 1 Beta 3 
[54]. We used Jeffreys-type priors for pi and theta, a bumin of 
100,000 sweeps and run 1,000,000 sweeps afterwards. Moreover, 
we calculated a Principal Component Analysis (PCA) using the 
adegenet package [55] for R 3.0.3 [56] replacing missing data by 



the mean frequency of the corresponding allele. For this analysis, 
we assigned all individuals to their most likely STRUCTURE 
clusters. Wild-caught individuals that clustered with domestic cats 
in STRUCTURE as well as captive individuals that clustered with 
wildcats in STRUCTURE were assigned to a unique cluster. We 
then simulated 100 first generation hybrids as well as 100 
backcrosses with either domestic cats or wildcats and performed 
a second PCA to inspect the degree of overlap of captive 
individuals with the simulated hybrids. 

For comparing genetic diversity, we calculated the mean 
number of alleles, allelic richness and inbreeding coefficients 
(Fis) for the captive, domestic and wild samples as well as for the 
genetic clusters in Fstat. GenAlEx was used to determine the 
expected (H^) and observed (HJ heterozygosity for each locus and 
each population and to test for deviations from Hardy-Weinberg 
ecjuilibrium (HWE). Furthermore, we tested the number and 
frequency of private alleles for domestic cats and wildcats 
(excluding potential hybrids) and analysed the number of alleles 
shared with captive cats as well as their mean frerjuencies. 

Effective population sizes (N,.) for each the captive population, 
domestic population and the wild population (excluding the 
hybrids) were assessed in ONeSAMP [57], which uses a Bayesian 
approach. The upper and lower bounds of the prior distribution 
for Ne were 2 and 500, respectively. The captive and wild samples 
were tested for genetic signatures of population bottlenecks using 
the Wilcoxon signed rank statistic implemented in Bottleneck 
1.2.02 [58]. We examined three different mutational models; the 
infinite alleles model (lAM), the stepwise mutational model (SMM) 
and the two-phase mutational model (TPM). However, we mainly 
considered the TPM, as this is the most likely mutation model for 
microsatellites [59]. 

The mean relatedness between individuals within the wild and 
captive population was calculated using Coancestry [60]. We 
assessed the performance of all seven relatedness estimators by 
comparing the results with known relationships and chose the 
estimator based on Wang [61], which showed the smallest 
deviation from known relationships and the smallest variance. 
Additionally, we tested the relatedness between the 19 known 
breeding pairs of the ex situ population. To measure individual F, 
four estimators are available in Coancestry. We used the TrioML 
estimator based on Wang [62] as it fitted our data best based on 
the relatedness data for the whole captive population (see above). 

Sequencing. DNA sequences were corrected and aligned by 
eye as no indels occurred. We excluded ambiguous data from the 
beginnings and ends of the fragments in the analyses. The final 
alignment contained 828 bp (positions 12642-13469 based upon 
GenBank entry FCU20753). Sequences were deposited in 
GenBank under the accession numbers KM246612-KM246625. 
The identification of haplotypes was carried out in DnaSP 5 [63] . 
For calculating a haplotype network we included 198 sequences 
(Table S4 & S5 in File SI). To verify the assignment of haplotypes 
to wildcats and domestic cats, we also performed a second analysis 
including reference samples for European wildcats and domestic 
cats published by DriscoU et al. [32] (identical haplotypes are given 
in Table S5 in File SI). The program TCS 1.18 [64] was used to 
construct a parsimony-based network representation of the 
mitochondrial haplotypes with a connection limit of 95%. 

Results 

Genetic structure 

Evidence for null alleles was detected only for single markers in 
single genetic clusters. All pairwise tests for linkage disequilibrium 
were non-significant (p>0.05). In the STRUCTURE analysis, we 
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found a clear distinction between individuals from the captive (red 
cluster in Fig. la) and wild (green cluster in Fig. la) populations at 
K = 2. AH domestic cats and all wild-caught domestic cats were 
assigned to the captive cluster (Fig. la). At K = 3 (optimal K value 
according to AK; Fig. S3 in File SI) the three groups "wild" (green 
cluster in Fig. lb), "captive" (blue cluster in Fig. lb) and 
"domestic" (red cluster in Fig. lb) were clearly distinguished. 
The most likely number of genetic clusters (K) revealed by 
Structure was six (Fig. Ic), reflecting geographical groups of 
breeding lines (Table 1). For example wildcats assigned to cluster 6 
(light blue in Fig. Ic) all originated from Scandinavian zoos or had 
ancestors from these, cluster 3 (yellow in Fig. Ic) comprised mainly 
ex situ wildcats stemming from one German zoo, and cluster 4 
(dark blue in Fig. Ic) mainly included individuals from eastern 
European zoos. The wild population (green in Fig. Ic) represented 
a single cluster together with six captive samples. These six 
samples came mostly from zoos within the range of the wild 
populations and one of these individuals was a foundling. However 
not all of them had wildcat mtDNA haplotypes. Another cluster 
(red in Fig. Ic) included all domestic samples together with the 
wild-caught domestic cats, some of which were morphologically 
suspicious. 

Gene-tic differentiation based upon Fg-p was significant for all 
pairs of genetic clusters (p = 0.001, Table 2), but not for Rst- The 
wild cluster (cluster 2) showed a higher differentiation to the 
captive clusters (mean Rst = 0.032, mean Fst = 0.187) than the 
domestic cluster (mean Rst = 0.015, mean Fst = 0.163). For the 
domestic cluster, the highest differentiation occurred to the wild 
cluster (Rst = 0.059; Fst = 0.254). The nested AMOVA based on 
FsT revealed that most of the genetic variation occurred within 
individuals (80%). Yet, a statistically significant portion was 
explained by the difference between the six genetic clusters 
(12%) and between the captive, domestic and wild samples (7%). 
In the RsT"based AMOVA, the population level (i.e. the clusters) 
explained no variation at aU, Imt the three categories (i.e. domestic, 
wild, captive) explained 7% of the variation. With Rst only 14% 
of the genetic variation occurred within individuals. The private 
alleles analysis revealed that captive cats shared 30 private alleles 
only with domestic cats (mean allele frequency of 2 1.0% ±3.1% 
SE) and six only with wildcats (8.6% ±4.2% SE). 

When analysing the data in NewHybrids, only two captive 
individuals were assigned as purebred wildcats, whereas most 
captive cats were either assigned as F2 hybrids (62 individuals) or 
backcrosses with domestic cats (13 individuals). The projected 
inertia of the PGA was 15.7% for the first axis and 6.9% for the 
second axis. The first axis separated the domestic cats (negatively 
loaded, cluster 1 in Fig. 2a) from the wildcats (positively loaded, 
cluster 2 in Fig. 2a). The captive clusters (3-6) had intermediate 
positions with negative loadings on the first £ixis. Captive 
individuals that clustered with the wildcats in the STRUCTURE 
analysis (cluster 7 in Fig. 2a) had a strong o\'erlap with the wildcat 
cluster (except for two individuals). Wild-caught indi\'iduals that 
had been assigned to the domestic/hybrid cluster in STRUC- 
TURE had a strong overlap with the domestic cluster in the PGA 
as well (cluster 8 in Fig. 2a). When performing the same analysis 
including the simulated hybrids and back-crosses, there was a 
substantial overlap of the captive clusters with first generation 
hybrids (clusters 3-6 in Fig. 2b) and of the wild-caught, 
morphologically suspicious individuals with second generation 
backcrosses with domestic cats (cluster 8 in Fig 2b). 

We found 14 mtDNA haplotypes in our samples (Fig. 3), two of 
which (WCl, WC2) were found in 31 in situ wildcats, 15 captive 
individuals as well as the European wildcat reference samples 
published by DriscoU et al. [32]. Further ten haplotypes (DGl- 



DGIO) were found in domestic cats (including the domestic cat 
reference samples), six known or suspected in situ hybrids and 45 
zoo samples (as well as GenBank sequences from domestic cats). 
One haplot)'pe found in four captive individuals (WC3) was 
identical to a haplotype found in an Iberian wildcat (Fsi257) by 
DriscoU et al. (2007). The remaining haplotype (WC4*) was found 
in 36 wild individuals from the Harz population as well as two 
captive individuals. It was closely related to the most common 
domestic haplotype DCl (two substitutions). Altogether 68% 
(n = 45) of the captive individuals had domestic haplotypes, 
whereas 32% (n = 21) had haplotypes found in wild populations 
of F. s. silvestris (Figs. Id, 3). Only one captive cluster 
(STRUCTURE cluster 6) was littie affected by mitochondrial 
introgression with only a single individual possessing a domestic 
haplotvpe. However, all individuals belonging to this cluster were 
assigned either as F2 hybrids or backcrosses with domestic cats in 
the NewHybrids analysis. 

Genetic diversity in captive and wild population 

The effective population sizes of the captive, wild and domestic 
populations were 70±1.74, 32±1.11 and 53±3.02 (mean ± SE), 
respectively. The captive clusters had a higher allelic richness and 
higher expected and observed heterozygosities than the wild 
cluster (Table 3). There was httie variation in genetic diversity 
among the captive clusters. 

The captive population showed significant evidence of a recent 
botdeneck under both lAM and TPM, but not under SMM. For 
the domestic population there was only significant evidence of a 
bottleneck under the lAM, whereas the wild population showed no 
signs of a botdeneck (Table 4). When the genetic clusters were 
analysed separately, two captive clusters showed no evidence of a 
bottleneck, whereas all other clusters showed evidence of a 
botdeneck under the lAM. 

The mean relatedness between dyads was r = — 0.044±0.007 
(mean ± SE) in the wild population and r = -0.042 ±0.004 in the 
captive population. The mean individual inbreeding coefficients 
were 0.149±0.019 in tiie wild and 0.139±0.019 in the captive 
population. In both populations most individuals (61%) had an 
individual inbreeding coefficient <0.125. The number of individ- 
uals with an inbreeding coefficient between 0.125 and 0.25 was 
18% and 14% for the wild and captive population respectively and 
the number of individuals with values above 0.25 was 23% and 
25%. 

Within the breeding pairs we found a mean r of 0.068 ±0.059. 
Most breeding pairs (47%) had a low relatedness r£0.05. 
However, there were five breeding pairs which showed relatedness 
values above the half-sib level (0.25) with two of them almost 
reaching the fuU sib level (0.5). 

Discussion 

Genetic structure of the European wildcat populations 

We found four distinct genetic clusters within the captive 
population of F. s. silvestris (Fig. Ic), which reflected their origin 
from certain zoos or regional groupings of zoos. The differenti- 
ation of the clusters is probably caused by the separation of 
breeding lines of groups of zoos that often exchanged individuals 
or zoos that were successful breeders. The existence of isolated 
breeding lines confirms a lack of knowledge in zoos concerning the 
stocks of other zoos. Therefore, an exchange of individuals is 
usually restricted to known partners. When testing the data in 
STRUCTURE, only six captive individuals were assigned to the 
wild cluster, some of which were foundlings from the in situ 
populations that had recentiy been integrated into the captive 
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Figure 1. Genetic clusters obtained from tKie Structure analysis (n = 189). Each individual is represented by a single vertical line, divided into 
K colours. The coloured segment shows the individual's estimated proportion of membership to that genetic cluster; (a) assignment for K = 2 (red = 
domestic/captive, green = wild); (b) assignment for K = 3 (red = domestic, green = wild, blue = captive); (c) assignment for K = 6 (cluster 1: 
domestic; cluster 2: wild, clusters 3-6: captive), (d) Analogue to the cluster assignments the haplotype affiliation for each individual is indicated by a 
vertical bar (orange = domestic haplotype, green = wildcat haplotype, white = missing data). For the definition of clusters (geographic origin and 
ancestry of individuals) see Table 1. 
doi:1 0.1 371 /journal.pone.01 06083.g001 



population. However, only four of them carried wild mtDNA 
haplotypes. At K = 2 we found a clear distinction between wild 
individuals on the one hand and captive individuals together with 
domestic cats on the other hand. When analysing the data in 
NewHybrids, only two captive individuals were assigned to 



purebred wildcats and in the PCA there was a strong overlap of 
the captive clusters with simulated hybrids and backcrosses. 
Furthermore, the level of dilferentiation based upon Rgf values 
was much lower between captive clusters and domestic cats than 
between captive clusters and wildcats. This suggests that the 



PLOS ONE I www.plosone.org 



5 



August 2014 | Volume 9 | Issue 8 | el 06083 



Genetic Integrity of Ex Situ Wildcats 



Table 1. Structure based clusters combined with available data on origin and ancestry of the included individuals. 





cluster 


grouping 


individuals lived in or originated from zoos in 


cluster 1 


domestic/ hybrid 


all domestic samples and some wild-caught (morphologically suspicious) individuals 


cluster 2 


wild 


wild-caught individuals from the Harz population and single ex situ samples from some German zoos 


cluster 3 


single German zoo 


almost exclusively individuals originating from one German zoo; all four individuals with the Iberian haplotype (WC3) 
were assigned to this cluster 


cluster 4 


east-european zoos 


cats originating from eastern European zoos or with ancestors from eastern European zoos 


cluster 5 


various zoos 


various zoos 


cluster 6 


skand. zoos 


cats originating from Scandinavian zoos or with ancestors from Scandinavian zoos 
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captive population is strongly influenced by gene flow with 
domestic cats and that the captive cats are not representative for 
the wild population. 

In the wild population we found two distinct mtDNA clades, 
which is in agreement with other genetic studies [65,66]. One 
mtDNA lineage (consisting of the haplotypes WCl, WC2 and 
WC3) also included wildcat samples from Rhineland-Palatinate 
and the Saarland as well as wildcat reference samples of DriscoU et 
al. [32] (see also table S5 in File SI). This lineage ("wildcat", 
Fig. 3) was well separated from the other clades and probably 
represents the autochthonous mitochondrial lineage of the 
European wildcat (including the substantially differentiated 
Iberian clade WC3). The second in situ lineage ("wildcat 
introgressed", WC4*, consisting exclusively of individuals from 
the Harz population and two captive individuals) was closely 
related to a common domestic cat mtDNA haplotype (DCl). We, 
therefore, suspect that this lineage originated by ancient 
introgression from domestic cats. As this mtDNA haplotype is 
already derived from the domestic cat haplotype and no domestic 
cat shared exactly the same haplotype, we suppose that 
introgression probably occurred hundreds to thousands of years 
ago (due to the relatively small genetic distances a molecular clock 
is not applicable). This scenario is a reasonable explanation as the 
wild population in the Harz region has had a long history of 
geographic isolation due to the extinction of most adjacent 
populations [26,65,66]. During this period of isolation the 
population has most likely declined due to persecution and some 
of the remaining wild individuals might well have mated with 
domestic cats. On the other hand, it remains unknown if this 
lineage is even more widespread in eastern or south-eastern 
Europe as samples from this region are scarce. Due to the 
differences to pure domestic cat haplotypes, we treat this 



introgressed in situ lineage as a second wildcat lineage, which is 
also representative for a natural population. 

In the wild population we found only few individuals (most of 
which were roadkHls that were known or suspected to be hybrids) 
who shared haplotypes with domestic cats, confirming that 
hybridization stiU occurs in the wild but at a low frequency. 
Hybridization with domestic cats has been detected in several free- 
ranging populations of the European wildcat [25,26,32,65-69]. In 
contrast to the past introgression event in the wild population, the 
captive population seems to be strongly affected by recent 
introgression, as two thirds of the individuals shared haplotypes 
with domestic cats (Fig. 1 and 3). We infer, based on the data 
presented here, that some founders of the captive population were 
hybrids. Zoos and wildlife sanctuaries are often confronted with 
stray kittens, which are suspected to be wildcats and are integrated 
into the captive population (as was the case with some of our 
sampled individuals). However, there is a high chance that many 
of these individuals are of hybrid origin. Humans are more likely 
to encounter litters or stray kittens of domestic cats, as these are 
tamer and have a higher probability of living close to human 
settlements [70]. Hybrids are difficult to identify based on their 
morphology and even the discrimination of some domestic cats 
from pure wildcats remains difficult [24,32,67,71]. Without 
genetic screening there remains a high risk of integrating hybrids 
or backcrosses into the captive population. 

Genetic integrity of the captive population 

Our results show that the captive population of the European 
wildcat represents a hybrid swarm. Only six captive individuals 
(7.7%) were assigned as wildcats in the STRUCTURE analysis (at 
K = 2), two of which had domestic cat mtDNA haplotypes, and 
only two individuals were clearly assigned in the analysis with 
NewHybrids (one with a domestic cat mtDNA haplotype). 



Table 2. Fsj values (above diagonal) and Rsj values (below diagonal) for genetic differentiation between the genetic clusters. 



Cluster 1 (domestic) Cluster2 (if/VaO Clusters Cluster 4 Clusters Clusters 



Cluster 1 




0.254* 


0.185* 


0.150* 


0.147* 


0.169* 


Cluster 2 


0.059* 




0.234* 


0.183* 


0.147* 


0.184* 


Cluster 3 


0.000 


0.033 




0.122* 


0.132* 


0.181* 


Cluster 4 


0.023 


0.019 


0.100* 




0.108* 


0.099* 


Cluster 5 


0.036* 


0.050* 


0.114* 


0.161* 




0.114* 


Cluster 6 


0.000 


0.024 


0.068* 


0.121* 


0.074* 





The in situ samples were assigned to cluster 2, whereas the domestic samples were assigned to cluster 1. Significant values are marked with * (P<0.05). 
doi:l 0.1 371 /journal.pone.Ol 06083.t002 
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Figure 2. Plot of the first two axes of tiie Principal component analyses (a) including all obtained genotypes, (b) including 
simulated first generation hybrids as well as backcrosses with either domestic cats or wildcats. Numbers represent the cluster numbers 
obtained in STRUCTURE; colours in (a) correspond to the clusters in Figure 1. Captive individuals that were assigned to the wildcat cluster were 
assigned to a unique group (7) as well as wild-caught individuals that were assigned to the domestic/hybrid cluster (8). The inserted graph shows the 
distribution of Eigenvalues. 
doi:1 0.1 371 /journal.pone.01 06083.g002 
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Altogether, 68% of the captive population possess domestic cat 
mtDNA haplotypes. Therefore, the captive population cannot be 
recommended for further breeding and reintroduction. As 
mtDNA is maternally inherited, the presence of domestic cat 
haplotypes in two thirds of the ex situ population does not allow 
any conclusions on the total amount of domestic cat nuclear 
introgression into the captive wildcat gene pool. Furthermore, the 
amount of introgression from male domestic cats remains 
unknown, as no Y-chromosome markers have been studied. 
However, the European wildcat is stUl widespread and it is possible 
to obtain new founders from the wild. Particularly the large 
populations in the Hunsrueck Mountains (Rhineland-Palatinate/ 
Saarland) currently show no indication of past introgression, but 
more detailed studies are required for any final conclusions. 
Individuals of hybrid origin cannot be recommended for breeding. 
This leaves many zoos with the difficult decision of what to do with 



the individuals which do not qualify as "pure" wildcats. One might 
question whether the introgression affecting the captive population 
occurred recendy or if it reflects another trace of an ancient 
introgression event. Hertwig et al. [65] suggested that "pure" 
wildcats without any traces of past introgression do not exist in 
Central Europe. They hypothesized that repeated introgression 
from domestic cats, followed by a spread and diversification of the 
wildcat populations might have played a major role in the 
establishment of European wildcat populations. However, the 
domestic cat haplotypes found in the captive population were only 
found in few wild individuals, all of which were declared to be 
morphologically suspicious and most probably were wUd-caught 
domestic cats. Although we did not sample the captive individuals 
ourselves, we occasionally came across some individuals, which 
were also morphologically suspect. This supports our hypothesis 
that recent hybrids were integrated into the captive population as 




Figure 3. Unrooted parsimony networit of mtDNA Kiaplotypes based on an 828 bp sequence of ND5. The sizes of the circles roughly 
reflect haplotype frequencies in our samples. 
doi:1 0.1 371/journal.pone.01 06083.g003 
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Table 3. Genetic diversity and inbreeding in captive, domestic and wild populations and the Structure based clusters. 



population/cluster 


n 


n> 




H„ 


He 


F,s 


cluster 1 (domestic) 


32 


8.2 


5.9 


0.731 


0.741 


0.005 


cluster 2 (wild) 


79 


4.7 


3.3 


0.539 


0.544 


-0.012 


cluster 3 (captive) 


9 


3.8 


3.8 


0.711 


0.603 


-0.170 


cluster 4 (captive) 


17 


5.4 


4.6 


0.600 


0.691 


0.108 


cluster 5 (captive) 


33 


6.3 


4.6 


0.709 


0.697 


-0.033 


cluster 6 (captive) 


12 


4.1 


3.9 


0.717 


0.664 


-0.138 


mean captive clusters 


71 


4.9 


5.7 


0.684 


0.649 


-0.040 


captive 


77 


7.7 


6.9 


0.681 


0.738 


0.085 


domestic 


32 


8.2 


8.1 


0.741 


0.741 


0.064 


wild 


73 


4.1 


3.7 


0.532 


0.533 


0.021 



n= sample size; na= number of alleles; Ar= allelic richness based on the lowest samples size (n = 9 for the cluster-based analysis; n = 32 for the three sampled 
populations); Ho= observed heterozygosity; He= expected heterozygosity; F|s= inbreeding coefficient. 
doi:l 0.1 371 /journal.pone.Ol 06083.t003 



founders. Based upon these conclusions we recommend to abstain 
from any further reintroductions of captive European wildcats 
from the existing captive stock. 

It is very likely that the reintroduced populations of the 
European wildcat in parts of Germany (e.g. Bavaria, Hesse) 
represent hybrid swarms as well as they stem from the captive 
population studied. Removing the reintroduced populations from 
the wild is not feasible as the administrative effort and costs for 
such a project would be far too high and it is very unlikely that all 
individuals could be caught and screened genetically. Moreover, it 
remains unknown to what extent the reintroduced wildcat 
populations originated from the captive population or immigrated 
from wild populations. A genetic survey of the reintroduced 
populations is therefore needed to assess their status. As the 
European wildcat currently expands its range (for example in 
Germany) [22,65,72], further reintroduction programmes are not 
needed and also not cost-efficient. Nevertheless, regional conser- 
vation administrations still aim at reintroduction in areas which 
are currently not colonized naturally. 

Recommendations for a future breeding strategy 

Outbreeding is a major threat for the success of ex situ 
conservation programmes [5,7,73]. Some principle decisions on 



the fate of the captive wildcat population have to be made by tiie 
holders and zoo associations. The frrst question concerns the 
priority of keeping the European wildcat in zoos. Space is a limited 
resource in zoos and it might be a wise decision to focus on 
breeding highly endangered species rather than non-threatened 
ones. However, zoos also have an educational function, which 
legitimates the display of native wild fauna. Therefore, if the 
strategic plan of zoo associations comprises a permanent 
husbandry of European wildcats, large parts of the current stock 
need to be replaced with pure wildcats. A future captive 
population should reflect the genetic structure and diversity found 
in the wild, particularly in populations which are little affected by 
hybridization with domestic cats. A sporadic integration of wild 
individuals into the captive breeding stock has been recommended 
by several authors [74,75]. In the case of the European wildcat, 
zoos are in the lucky position that they can obtain individuals from 
healthy wild populations and, therefore, are able to increase the 
genetic variability of their breeding stock. Nevertheless, each new 
individual should be screened genetically before it is integrated 
into the population. 

The current plan to establish a studbook for the European 
wildcat also requires a coordinated breeding strategy. The most 
widely accepted management strategy is minimizing kinship [76], 



Table 4. One-tailed p-values of the Wilcoxon test for heterozygosity excess under three mutational models for captive population, 
wild population and the genetic clusters. 



population/cluster lAM TPM SMM 



captive 


0.00049* 


0.00244* 


0.98389 


domestic 


0.00488* 


0.53906 


0.98779 


wild 


0.04199 


0.31250 


0.98779 


clusterl 


0.00684* 


0.57715 


0.99316 


cluster2 


0.00684* 


0.24609 


0.90332 


clusters 


0.06543 


0.24609 


0.57715 


cluster4 


0.06543 


0.34766 


0.91992 


clusters 


0.00146* 


0.27832 


0.98389 


cluster6 


0.00244* 


0.00488* 


0.11621 



lAM: infinite alleles model; TPM: two-phase mutational model; SMM: stepwise mutational model. Significant values are marked with * (P<0.005). 
doi:l 0.1 371 /journal.pone.Ol 06083.t004 
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which aims at minimizing the overall level of relationship in the 
population and maximizing the retention of genetic diversity 
[77,78]. However, for this approach, detailed pedigree data is 
needed to find optimal breeding pairs [77,79]. Genetic analysis 
during the establishment of a new captive population can help to 
minimize the risk of inbreeding and outbreeding. An artificial 
fragmentation of captive populations into several more or less 
independent subpopulations has been recommended by several 
authors to retain a maximum of genetic diversity [74,80]. 
Although such a fragmentation has been achieved in the current 
ex situ population of F. s. silveslris just by chance (due to the 
missing network of holders), and might have possibly helped to 
maintain the strong genetic diversity, it is likely that this diversity 
has mainly been caused by admixture with domestic cats. 
Increased genetic diversity caused by admixture of several 
evolutionary hneages has also been reported from invasive non- 
native populations [35]. 

Genetic monitoring of ex situ conservation projects 

The importance of genetic assessments of captive populations 
has already been recognized decades ago [81]. Molecular gcn(^tic 
studies of captive populations have provided valuable information 
for the improvement of conservation breeding programmes [7]. 
These include the detection of hybrids [82], the detection of 
genetically valuable individuals that had been excluded from 
breeding [83], the vahdation of studbook data and detection of 
errors [11,12], and finally the detection of insufficient represen- 
tation of genetic diversity found in the wild [10]. 

The present study shows that genetic analyses provide valuable 
insight into the status of wild and captive populations. Using two 
different marker systems (microsateUites, mtDNA) helps to unravel 
introgression patterns. Future studies should also integrate Y- 
chromosome markers to detect potential male introgression. The 
integration of hybrids into ex situ populations can have drastic 
effects on the genetic composition of captive populations. 
Hybridisation with closely related species or domesticated 
congeners represents a severe threat to the conservation of 
endangered species both in situ and ex situ [35,82,84—87]. Until 
now, many individuals from the captive population of the 
European wildcat have been reintroduced to the wild without 
any test of their genetic integrity. Such reintroductions may 
strongly deteriorate the genetic integrity of wild populations and 
promote cryptic extinctions [85] . In some parts of Europe (e.g. 
Scotland) the wild populations of the wildcat are severely 
threatened by hybridisation [22]. Therefore, measures to assess 
and constrict any further admixture of wild and domestic cats are 
of great importance [22]. Although introgression must have 
occurred also in nature (perhaps even repeatedly), the European 
wildcat can still be regarded as a distinct taxon that exists in 
sympatry with the domestic cat and needs to be managed as a 
separate unit. 
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